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Recently, Glass and Krisch have extended the Vaidya radiating metric to 
include both a radiation fluid and a string fluid [1999 Class. Quantum Grav. 
16, 1175]. Mass diffusion in the extended Schwarzschild atmosphere was stud- 
ied. The continuous Solutions of classical diffusive transport are believed to 
describe the envelope of underlying fractal behavior. In this work we examine 



the classical picture at scales on which fractal behavior might be evident. 

O 

\o 

O ' PACS numbers: 04.20.Jb, 5.45Df, 64.6 

O ' 
O 

"o 
er- 

u 

bJQ. 

> 

• rH 

IS 



*Permanent address: Physics Department, University of Windsor, Ontario N9B 3P4, Canada 



1 



I. INTRODUCTION 



Since its discovery in 1916, the Schwarzschild vacuum Solution has been a source of many 
insights into relativistic physics. 

g s J ) h dx a dx h = (1 - 2m /r)du 2 + Idudr - r 2 dü 2 . (1) 

This Solution and its extensions contain new views of relativistic physics on a variety 
of scales, from the size of the Universe down to the size of a stellar black hole. When 
contemplating quantum effects at classical boundaries such as the Schwarzschild horizon, 
the scales become much smaller down to the Planck length L p = ^JhG/c 3 = 1.6 x 10~ 33 cm. 
The classical point particle electron may perhaps be resolved into Planck scale strings bits 
with size proportional to the Planck length and with 10 20 strings bits Atting into a classical 
electron radius. 

Vaidya U extended the Solution by allowing the Schwarzschild mass parameter to be 
a function of retarded time, m(u), and discovered a spherically Symmetrie null fluid atmo- 
sphere beyond the horizon at 2mo- Glass and Krisch (GK) [g] pointed out that allowing the 
Schwarzschild mass parameter to be a function of radial position and retarded time, m(u, r), 
creates an atmosphere with both a null fluid and string fluid stress-energy content described 
by metric g^ h '. The density of the string fluid, p, obeys a classical diffusion equation in the 
Schwarzschild spacetime: 

d u p = r~ 2 d r (Dr 2 d r p) (2) 

where D(r) is the variable diffusivity || found in anomalous diffusion. One may show [|J, 
if the mass is diffusing in a flowspace with ^J—g = r~^/(i?) (in ||] ß — 1 — S), that the 
functional dependence of the diffusivity is determined. 

d u m = r /3 9 r (Dr~ /3 9 r m) (3) 
D(r) = D r ß ~ 2 . 
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The continuous Solutions of classical diffusive transport are believed to describe the envelope 
of underlying fractal behavior ||. Fractal behavior is usually associated with mass scaling 
but the radial dependence of the diffusivity || carries additional information about the 
Connectivity of the underlying fractal Substrate, and so we wish to examine the classical 
picture at scales on which fractal behavior might be evident. These smaller scales are still 
classical (i.e. larger than quantum scales). We consider a 2-parameter family of spacetimes 
with metric function A— 1 — 2m(u,r)/r 

g s a t ale dx a dx b = Adu 2 + 2dudr - X 2 - a r a dQ 2 , (4) 

where A has units of length, and is required for spheres to have area of (length) 2 . This 
metric is capable of describing different scales for a particular central source. For example, 
a particular mass function m(u,r) with at = 2, could describe the Schwarzschild vacuum 
outside the source. For a ^ 2 it would describe the source at different scales where we 
expect evidence of fractal behavior to appear. Fractal dimensions are generally less than the 
Euclidean dimension, and we will describe an averaging process on 2-spheres for a < 2. The 
scale parameter a cannot be removed by a coordinate transformation. If one attempts the 
appearance of the usual area coordinate by y 2 = X 2 ~ a r a , then a will emerge in the (u,r) part 
of 

gsca/e_ M oreover; the manifestly coordinate-independent sectional curvature of the <p) 
two-surfaces, given in Eq.(|i~5|) below, provides a physically meaningful mass expression in 
terms of at. 

The interpretation of a as a fractal scale follows from considering diffusive processes in 
the atmosphere. If V(u, r) is a quantity that is diffusing, the diffusion equation is 

d u V = r- a d r {Dr a d r V) (5) 

and a + 1 is identified as the dimension of the fractal Substrate on which the diffusion is 
occurring ||, J7J. At smaller quantum levels, V can be a probability measure on the space 
of radial paths [|J. 

The fluid atmosphere described by this metric is discussed in the next section. In section 
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III diffusion and scale sizes are discussed. Definitions of topological and Hausdorff dimen- 
sions are given. Applications of the fractal view are presented in section IV, along with 
specific mass Solutions. The paper is summarized in the final section. 

In this work Latin indices ränge over (0,1,2,3) = (u, r, (p). Our sign Conventions are 
2^4 c; [afe] = A e R e cab , and R ab = R c abc - Overdots abbreviate d/du, and primes abbreviate 
d/dr. Overhead carets denote unit vectors. We use units where G = c = 1. Einstein's field 
equations are G ab = —8iiT ab , and the metric signature is (+,-,-,-). 



II. THE ATMOSPHERE 

Metric g s ^ le describes an atmosphere that admits a two-fluid description of matter with 
outward flowing short-wavelength photons (sometimes called a "null fluid"). 

With the use of a Newman- Penrose null tetrad the Einstein tensor for the atmosphere is 
computed from (|4]) and given by 

G ab = -2$ n (l a n b + nj b + m a rh b + rh a m b ) (6) 
-2% n a n b - 2<£>2 2 Uf> - GAg ab . 

Here the null tetrad components of the Ricci tensor are 



a(2 - a) 



_2m'-rm" a 2 \ a ~ 2 (a-2)(a + 2)m . . 

11 _ 4 r 2 lQ r 2 ^ 4 r a ^ g r 3 ' ' ' 



am n (2 - n ) . , 

* 22 = _ 2? + lft^ (1 " 2m/r) ' (7C) 



A = Jt/a4= rm» + 2(a-lK + («-2)(8a-2)m 
12r 2 24r 3 
a(3a -4) \ a ~ 2 



48r 2 ' 12r c 

The only non-zero component of the Weyl tensor is 

(a + 1)?77. (2 + a)m' — rm" at A Q_2 , , 

^2 = -- — — h - — — 1 • (8) 

3r 3 g r 2 12r 2 Q r a V ) 

The metric is Petrov type D with l a and n a principal null geodesic vectors 



l a dx a = du, l a d a = d r (9a) 

n a dx a = {A/2)du + dr, n a d a = d u - {A/2)d r (9b) 

m a dx a = -(l/ v /2)A 1 - a/2 r ö/2 (^ + i wy&dtp), m a d a = (l/v^A / 2 " 1 ^ / 2 ^ + —d v ). (9c) 

When a = 2 all quantities have the spacetime values for metric g^ b K discussed in |2j . 

In order to clearly see the matter content we introduce a time-like unit velocity vector 
v a and three unit space-like vectors f a , § a , (p a such that 

alT^ = VaVb ~ f a f b - d a d b - (f a (fi b . 

The unit vectors are defined by 

v a dx a = A 1/2 du + A~ 1/2 dr, v a d a = A~ 1/2 d u , (10a) 

r a dxa = A'^dr, f a d a = A~ 1/2 d u - A 1/2 d r , (10b) 

$ a dx a = \ x - at2 r a l 2 dd, td a = -\ a/2 - l r- a/2 d ß , (10c) 

a dx a = \ 1 ~ a/2 r a/2 smi9dip, <p a d a = -\ a/2 ~ 1 (r a/2 sm$)~ 1 d ip . (lOd) 

The Einstein tensor (|6]) yields the energy-momentum tensor 

T a b = IplJb + PVaVb + PrTah + P±(^a^b + fiafib) + q a Vb + V a q b (11) 



where 



am a (2 — a)A 2 

am' X a ~ 2 a(a — 1) (a — 2) r , , , . 

= — + — -^T 1 + - ~ 2m)], (12b) 



am' X a 2 a(a — 1) (a — 2) r 



nr 



" (a-2)m' (a-2) 2 m a(a - 2) 



8 ^ = " J ~ 2r 3 + -V^' (12d) 



r . a(2 — a), . / * 

87rg a = [A \ 2 ; f a . 12e 

The string equation of state p + p r = holds for the scale with a = 2. The decomposition of 
the energy-momentum tensor given in Eq. flllD is less compelling ona^2 scales since there 



is heat flow in addition to the null fluid. The heat flow arises from the $oo^a^b term in G a b. 
This is seen only on a ^ 2 scales and is the intermediate electromagnetic field part of the 
null fluid, which has only a far field part in the usual Vaidya large scale view. A different 
covariant decomposition using l a dx a = A~ l {y a — f a )dx a provides a Single fluid with radial 
heat flow. The energy-momentum tensor decomposes as 

Tab = PV a Vb + Prfah + P±('$a$b + fiafib) + q a V b + V a q b (13) 

where 

_ x am am! X a ~ 2 a(a — 1) (a — 2) 



8vrp = 


—A 


_ x am 


+ 


am 


8irp r = 


—A 


_ x am 




am' 


8vrg a = 


(A- 


l am. 

r 2 > 







-[Aam — a(r — 2m)], (14a) 



+ ^-T^ - hr^ [4öm + a{r ~ 2m)] ' (Mb) 

(14c) 



We note that 



_ .1 2am a(a — 2) 
MP+Pr) = -A l — ^(r-2m). 



In the static case the string equation of state p + p r = holds on the horizon for all a. For 
mass m(r) which falls off at least as fast as 0(l/r 2 ), the string equation of state also holds 
at future null infinity. 

The metric function m(u,r) measures the System mass only for a — 2, the GK metric. 
In general, spherical symmetry allows the mass within two-surfaces of constant u and r to 
be invariantly defined by the sectional curvature || of those surfaces: 

-2M/r 3 : = R abcd r<p h $ c d (15) 
a 2 m a 2 X a ~ 2 



+ 



When ol = 2 then M = m. 



III. DIFFUSION AND SCALE 
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A. Diffusion 



Continuous diffusive processes can envelope underlying fractal behavior. In the a = 2 
Schwarzschild atmosphere, the density p and the mass function m diffuse according to Eqs.(^) 
and ([|). By analogy with classical diffusion of density in the Schwarzschild spacetime, we 
assume that the density enveloping fractal behavior obeys 

d u T = r- a d r (Dr a d r T), (16) 

where the mass density T is defined on u = const surfaces, using the sectional curvature 
mass, by the relation M = J T^J—gdßx so that 

47r(r - r ) := X a - 2 r' a d r M. (17) 

The relationship between the energy density p, seen by observer v a , in (|1 2b| ) and the mass 
density is then given by 

4vrr = 4vrp + (— ) [(A/r) a ~ V(2m' - 1) - 2a(2m' - a + 1)] (18) 
+ (^)(A/rr 2 [(3 - a)(A/r)- 2 - 1] - ^^[8am - a(r - 2m)}. 

In the region where a = 2, the two densities T and p are equal. 



Assumption (|16D is equivalent to assuming a Fickian mass current of the form 

d u M = A-nD(r)r a d r T, (19) 

and then d r M = 47rA 2 ~ Q r a (r — T ) yields (|TBp. Furthermore, if the sectional curvature mass 
M diffuses as 

d u M = r ß d r (Dr- ß d r M) (20) 
then D(r) is determined by d u M in ( p!9]) and d r M: 

D(r) = D r ß ~ a . (21) 
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Note that if a = 2 the M diffusion equation is just Eq.(^) with ß = 1 — 8. For any a, when 
D = D , it follows that ß — at and that M and T diffuse reciprocally. 
Equation (PÖ| ) admits the homogeneous Solution 

M hom (r) = M + M 1 (r/X) 1+a (22) 

which can be added to each time-dependent Solution. There are two important observations 
that follow from these results. First, the metric of Mh om is not asymptotically ffat and the 
diffusive Solution would need to be matched to a Solution with appropriate asymptotic be- 
havior. Second, the diffusivity depends on a coordinate and not on the size of any correlated 
mass grouping. This implies that there may be multiple distance scales in the atmosphere. 



B. Dimensions 



To better understand fractal scales we brieffy discuss fractal dimensions ]nj • We cannot 
use the intuitive notion of the number of coordinates needed to locate a point in some region 
as the dimension of the region, since Peano constructed a curve (1 coordinate) which fills a 
Square (2 coordinates). Consider a point in set Q, with Q a subset of an Euclidean space 
of integer dimension dß- The notion of topological dimension involves the idea that neigh- 
borhoods around a point have boundaries of smaller dimension. Precisely, the topological 
dimension d^(fi) of set Q is the smallest integer such that each point has small neighborhoods 
with boundaries of lesser dimension. Since c?t(0) = — 1 it is clear that dy(fi) < dß- 



The Hausdorff dimension 11] of self-similar sets (fractals) is given in terms of coverings 
by open sets and is defined as 

d H {ü) := lim^^ 
r^o In r 

where N(r) is the smallest number of open balls needed to cover Q. The defmition of du 
is equivalent to the power law N(r) ~ const x r~ diI . When du = 3, N(r) corresponds to a 
3-dimensional Newtonian specific density. 
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The Hausdorff dimension of a fractal set is not always an integer. The scaling dimension 
ds of fractal sets is more intuitive and, since geometric fractals display regulär self-similarity, 
is equal to the Hausdorff dimension: 

d (Vi) ^ ( num ^ er °f re P ncas ) 
log(magnification factor) 

Two famous examples are ^(Cantor set) = in dß = 1, and <is(Sierpinski gasket) = 
|°^| in dß = 2. An integer valued example is <i,s(Sierpinski gasket) = in dß = 3, since 
all Sierpinski gaskets obey d s (Vl) = 1 ° g 1 ogf 2 | 1 ' > • ^ n g enera l> dr(£l) < d H (ü) < d E - 



C. Scales 

Maximum scale size is independent of any particular application. Since a + 1 is the 
dimension of the fractal Substrate, bounded by dß, « max = 2. 

Consider a vacuum Solution at a = 2 and ask how it is produced. Beginning at a small 
scale, one considers the atmosphere to be made of many regions. Within each region an 
average energy density is found; in some of the regions it will be negative, in some positive, 
and in some zero. Then a larger region, containing several small regions is considered and 
again an average energy density is produced until in the end we have the vacuum atmosphere 
with all of the fluctuations averaging to zero. This process is just a renormalization group 
averaging procedure and, just as in renormalization group calculations on a lattice, it is 
necessary at each averaging step to rescale. In lattice calculations, one rescales to maintain 
the lattice spacing. In an atmospheric calculation with a horizon, the rescaling is done to 
always keep the physical observers outside the horizon. There are several ways in which this 
could be accomplished. 

One method is to have an observer at constant r examine an area of his 2-sphere at 
different scales while maintaining the same coordinate value of r at all scales. The entire 
2-sphere at a = 2 can be tiled by elements of area 5A 2 . The sum of all the tiles is the 
total area of the sphere Airr 2 . At small scales, the observer chooses a region of the sphere 
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containing many area elements to average over. The result of this averaging must be rescaled 
at each step, as in lattice renormalization, so that the rescaled area is again 5A 2 . With this 
averaging procedure, an observer outside the horizon for a = 2 will remain outside of the 
horizon for all a. The a = 2 horizon position will persist at all scales. 

To see the effect of the rescaling, consider a typical averaging process. Assume that a = 2 
is the vacuum Solution and let the observer look at an area 5A 2 and solid angle 6Q 2 with 
5A 2 = r 2 5Vl 2 . This region is the result of averaging over larger regions at the lower scale and 
then rescaling the area size. The larger region would have an area 5A a = \ 2 ~ a r a 5Vl a . The 
rescaling would adjust 5Q a to the size of 5Q 2 and thereby rescale the areas. The rescaling 
formulas are 

(\/rf- a 5VL a i — ► 5Q 2 (23) 
5A a i — > 5A 2 . 

The scaling factor (A/r) 2 ~ ö is expected to be < 1. For metrics without a horizon, the scale 
factor can be chosen. When there is a horizon, the scale factor is determined by its value 
on the horizon. To explore the numerical values of this factor in metrics with a horizon, a 
specific mass Solution is needed. 

IV. APPLICATIONS 

A. Constant mass Solution 

A specific mass Solution is necessary in order to explicitly treat diffusion and fractal 
behavior. We have seen above that the sectional curvature mass can be written as 

2 

2M =^(2m-r)+\ a - 2 r 3 - a . (24) 

To illustrate the rescaling process, we explore the case of constant M. For a = 2 there is a 
Schwarzschild vacuum with an horizon at r = 2m . Assuming that all scales see the horizon 
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at 2mo we can determine the scaling factor. From Eq.fl24|) with M = mo, and r = 2m = 2mo 
it follows that 



(2m /A) 



OL-2 



or A = 2mo- 

With m(r) from Eq.(j24j) the fluid parameters in Eq.flUD are 



8vrg a = 0, (25a) 

6(a-2) / . 

87rp= — -(m r- 3 ), (25b) 

et 

8Tr p r = - 2(a ~ 2) [m r- 3 + A- 2 (A/r) Q ], (25c) 
et 

8ttp ± = (— ) 2 [-2m r- 3 + A~ 2 (A/r) Q ]. (25d) 
a 



For this simple example, p, resulting from fluetuations, is negative [I2| and there is no 
bounding p r = surface. The radial profile describes an atmosphere asymptotically ap- 
proaching vaeuum. If one examines the at < 2 metric, it appears not to be asymptotically 
flat. However, asymptotic flatness applies only to the observer's a = 2 spacetime. 



B. Solution with finite boundaries 

The M = mo Solution is part of the homogeneous Solution, Eq. (|22|) . The complete 
homogeneous Solution, modified to produce vaeuum Schwarzschild with a horizon at r = 2m 
is 

M = m + m 1 (2-a)(r/X) a+1 (26) 

with 

4m 4m! (2 — a) , ,, N ^,i r 2A . „ 
a 2 a, 

Evaluating this equation on the horizon at r = 2mo to find the scaling factors gives the 
equation, with x = X/2m , 
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1 + (7rn/mo)(2 - a)x- {a+1) - x a - 2 = 0. (28) 

This more complex scaling equation has no Solution for ra\ > m , and may describe Solutions 
without an horizon. For ra\ « rriQ there are two scaling factors. Some typical values for 
mi/mo = 0.01 are 

a = 7/4, x ~ 1, 1/7.6 (29) 
a = 6/4, x~ 1, 1/11.9 
a = 5/4, x ~ 1, 1/24.5 
a = 1, x ~ 1, 1/99 

For the complete homogeneous Solution the fluid parameters are 

8nq a = 0, (30a) 
Svrp = 2(a ~ 2) [3m r" 3 + (4 - 5a) (mx/A 3 ) (r/A) a ~ 2 ] , (30b) 

OL 

8vrp r = - 2(Q ~ 2) [m r- 3 + A~ 2 (A/r) a - 3a(m 1 /A 3 )(r/A) Q ~ 2 ], (30c) 

OL 

87r P± = (— ) 2 [-2m r 3 + A- 2 (A/r) Q ] + (^^)[(5a - 4)(2m!/r 3 )(r/A) Q+1 ]. (30d) 

OL OL 

For Solution ([27]) there are boundaries for some a. The boundary radii follow from p r = 
at A = 2m . With y = rb/2mo, 

1 + 2y 3 - Q - 3a(m 1 /m )y a+1 = 0. (31) 

Again choosing mi/m = 0.01, there are boundaries which get larger as a decreases: 

a = 7/4, y ~ 11.5 (32) 

a = 6/4, y ~ 45 

a = 5/4, y ~ 2800 

a — 1, no boundary. 

Since a = 2 is a vacuum Solution, one could attempt construction of a matching vacuum for 
the lower values of a. However, the lack of a boundary at a — 1 implies that the boundaries 
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mark a change from pressure to tension rather than a boundary to vacuum. This supports 
a picture of an atmosphere which extends farther as the size of the observed fluctuation 
decreases, all asymptotically approaching vacuum. We again remark that asymptotic flatness 
applies only to the observer's a = 2 spacetime. 

C. THE SIERPINSKI GASKET 

These ideas can be applied to diffusive processes whose continuum Solutions are expected 
to envelope fractal behavior at smaller scales. As an example, we will use some existing 
discussions of enveloped fractal structures in flat Spaces to discuss the approximate sizes of 
ol and ß where D(r) = D r^~ a = D r @ . The numerical parameters of interest are a, the 
diffusivity factor 0, the Euclidean embedding dimension dp in which the fractal Substrate 
is embedded, and a fractal dimension dp (which may not be djj). Alexander and Orbach || 
and O'Shaughnessy and Procaccia || provide some numerical values of dp and for various 
fractal processes. The fractal dimension dp is defined through the diffusion equation which 
is written as 

dY 1 öpor-V^örr) 
du r d F -^ Q r 

In dp = 3, the fractal dimension for the Sierpinski gasket is dp = 2 and = 0.59. The 
corresponding values of a and ß are a — 1, ß — 1.585. For dp = 4, the fractal dimension is 
dp = 2.31 and = 0.81. We find a = 1.32 and ß = 2.13. 

For all of the above examples, we could find a u-independent Solution of the diffusion 
equations as 

r = T r^, 
M = .U„ • .!/;/•" '. 

The mass scaling is set entirely by the metric parameter a, while the density depends on 
the scaling of the diffusivity. 
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V. CONCLUSION 



Introducing a parameter, a, into the extended Schwarzschild metric, produces an at- 
mosphere with processes that can be linked to a fractal mass distribution. Specific values 
of a can be associated with particular fractal Substrates in a small scale region of the at- 
mosphere. The continuous Solutions of an anomalous diffusion equation in the atmosphere 
envelope smaller scale fractal processes. Anomalous diffusion can sample fractal processes 
on a topological level, giving Information about the Connectivity as well as the fractal di- 
mension. Fluid parameters from the field equations can be studied for a ränge of a, allowing 
the effects of fluctuations from one scale to another to be examined. 

Other models are possible. For example, a fractal Schwarzschild cell could be used in a 
universe model like the Lindquist-Wheeler model | 13fl , perhaps describing the fractal galactic 
density distribution. If the galactic distribution can be modeled by a Single fractal process 
in which both mass and density diffuse, then on time scales that are short compared with 
the diffusion time, the mass would scale with the fractal dimension while the scaling of the 
density would involve topological elements of the fractal Substrate. In this paper only a 
single level of scale reduction has been studied. It is possible that an ensemble of Solutions, 
a multifractal set summed over a, could be used to model more complex physical Systems. 

Another direction for future work is a possible link between fractal atmospheric structure 
and chaos. A variety of perturbations to a black hole Solution have been shown to lead to 



chaotic behavior, most often chaotic particle orbits [|14[]. For example, Bombelli and Calzetta 



15fl have shown that gravitational perturbations of the Schwarzschild Solution can produce 



chaotic orbits. In general, the effects of chaos ränge from unusual production of gravitational 



waves |16[ to a possible brightening of the black hole |17| . Orbiting particles with spin may 



also exhibit chaotic behavior []i~8 |, |fL9 |. A relation between a phase space fractal dimension 
and chaos has been established for multi black hole spacetimes [^ÜJ and for black hole chaotic 



scattering processes piß , f22|| . We have considered an atmosphere that has fractal structure 
in the physical spacetime. At one scale, the Solution can describe an ideal black hole, or 
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with time dependence and space dependence, a black hole with a two fluid atmosphere. The 
effect of perturbations of this atmosphere and the effects of the physical fractal structure on 
particle orbits is an interesting area for investigation. 
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